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The existing moving average control charts can be only applied when all observations in the data are 
determined, precise, and certain. But, in practice, the data from the weather monitoring is not exact 
and express in the interval. In this situation, the available monitoring plans cannot be applied for 

the monitoring of weather data. A new moving average control chart for the normal distribution is 
offered under the neutrosophic statistics. The parameters of the offered chart are determined through 
simulation under neutrosophic statistics. The comparison study shows the superiority of the proposed 
chart over the moving average control chart under classical statistics. A real example from the weather 
is chosen to present the implementation of the chart. From the simulation study and real data, the 
proposed chart is found to be effective to be applied for temperature monitoring than the existing 
control chart. 


The monitoring of the process is done to reduce the percentage of non-conforming items during the process. 
The control charts are very important tools and used to watch the effect of extraneous factors on the process. 
The timely alert about the shift in the process is helpful to identify the reasons for this shift. The Shewhart 
control charts are simple to apply in the industry but detect only the larger shift in the process. Therefore, the 
moving average (MA), cumulative sum (CUSUM), and exponentially weighted moving average (EWMA) are 
the alternatives of the Shewhart control chart. These charts are sensitive to detect a small shift in the process. A 
little attention has been paid to design MA control charts for various situations. Chen and Yang! developed an 
economic model for the MA chart. Wong et al.* discussed the sensitivity power of the MA chart. Khoo and Wong? 
and Areepong* proposed double MA control charts. Mohsin et al.? worked on a weighted MA chart using loss 
function. Alghamdi et al.° designed the MA chart for the Weibull distribution. Ye et al.’ applied quality control 
methods for air data. Su et al.? used the machine learning technique for metrology data. 

The control charts under classical statistics can be only applied when all observations from the industrial 
process are determined. In complex systems such as monitoring the level of water in a river, the process data may 
not be determined. “Fuzzy data exist ubiquitously in the modern manufacturing process”’. Therefore, an alterna- 
tive to the existing chart is the use of a control chart using the fuzzy approach. A piece of detailed information 
about the fuzzy control charts can be seen in references!? 76, 

[he fuzzy logic that provides information about the measure of truth and falseness is the special case of 
neutrosophic logic, see Smarandache'”. The neutrosophic logic considers an additional measure is known 
as the measure of indeterminacy. A detailed discussion on neutrosophic logic can be seen in references! ?*. 
Smarandache? introduced the neutrosophic statistics (NSO using the idea of neutrosophic logic. The NS is a 
generalization of the classical statistics and deals with the data having Neutrosophy. Chen et al.°°°! introduced 
the neutrosophic numbers with applications. The np and X-bar control chart using NS were developed by Aslam 
et al.°?°?. Aslam et al.22% worked on EWMA charts under NS. 

A rich literature is available on control charts under classical statistics and fuzzy approach. By exploring the 
literature and best of our knowledge, there is no work on moving average control chart under NS. In this paper, 
we will introduce this neutrosophic moving average (NMA) chart and evaluate its performance over the existing 
MA charts under classical statistics. We will present a real example from the weather-monitoring department. 
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Materials and Methods 

Suppose that Xy = Xz + bwIw; Ixe[Ir, Iy]be a neutrosophic random variable consists of the variable under clas- 
sical statistics Xr and indeterminate part by Iw; In€[Iz, Iv]. The neutrosophic Xye[Xr, Xy |; Ine [Iz Iv] reduces to 
X when I, =0. Suppose ny« [nr, ny] presents the neutrosophic group size. Suppose that X ¡ye IX iX iu | denotes 
the neutrosophic sample average for ith subgroup. Here we assume that Xj follows the neutrosophic normal 
distribution with a population mean uye[uz, M1] and variance ope lof. on for? = ,,2,...,and J = 1,2,.+<57- 
By following Montgomery”? NMA statistic is defined as 


WN 


MAin = ; wNe[wr, wu] (1) 

where wwe [wr, wy] shows the span at a time i. Note here that the statistic is given in Eq. (1) is similar to the 

exponentially weighted moving average (EWMA) statistic. Ihe main difference between MA statistic and the 

EWMA statistic is the sensitivity that each statistic shows in the calculation of the data. Ihe EWMA gives the 

higher weights to the current values while the MA gives equal weight to all values in the data, see Li et al.7. 
The MAjne[MAjz, MAjy] statistic in neutrosophic form can be written as 


MA¡n = MA + cn IN; IvelIz, Iu] (2) 


where MA is the determined part and cyIy; Inel, Iy] is indeterminate parts of MAjne[MAjz, MA¡y]. The 
MA statistic mentioned by Montgomery” is a special case of the proposed MA¡ye[MA¡L, MA;¡y]. The proposed 
MAiN€[MAir, MA¡y ] becomes MA statistic if I, =0. The neutrosophic mean and variance of MA¡ye[MA;¡¡, MAiu] 
when the process is an in-controli > wy; wye[wL, wy ] are given as follows 


EN (MAiN) = Mon; MA¡ye[MA¡L, MAiu]. one [Mol Hou] (3) 


and 





2 
pr ON . 2 2 2 
VN(MAiN) = re ; nue [ni nu], we[wL, wu]; owe |o. og] (4) 
N 


Based on the given information, the proposed control chart is mentioned as below. 
Step 1: Select a neutrosophic sample of size nwe[nr,nu] from the production and compute 
MA¡ne[MA¡L, MAiu |for specified wwe [wr, wy]. 


Step 2: Declare the process in-control if LCLy < MAjn < UCLy. 
The operational process of the proposed control chart depends on two neutrosophic control limits, which 


are given as 
ON 
LCLy = Mon — kw 4] — —5 Wye€[wL, wu], nyelni nu] (5) 
HNWN 


UCLw = Hon + kn ; WN€[wr, wu], ny e[nr, ny] (6) 





AN WN 


where ky e [kr, ky] is neutrosophic control limits. 


Neutrosophic Monte Carlo simulation 

In this section, we introduce the neutrosophic Monte Carlo (NMC) simulation for the proposed NMA control 
chart. Suppose that be neutrosophic mean for in-control process and Win = Mon + consu in€[ ui, Miu |, where 
cis a shift constant. Let row e [ror roy] be the specified neutrosophic average run length (NARL) when the process 
is in-control state, for more details, the reader may refer to Li et al.%. The NMC simulation is stated as follows. 


Step-1: Generate a random sample of size ny«e [nr, ny] from the neutrosophic standard normal distribution 
with Lov € [tor., ou] and variance ofe [of, 0 |. Compute Xine |Xir, Xiu | for ith subgroup. 

Step-2: Compute the statistic MAjne[MAjz, MA¡y] and plot it on LCLye[LCL¡, LCLy] and 
UCLye[UCL¡, UCLy]. Note the first out-of-control value, which is called the run length. 

Step-3: Repeat the process 10,000 times and compute the neutrosophic average run length (NARL) and 
neutrosophic standard deviation (NSD). Choose kye[kr, ky] for which NARL for in control process, say 
ARLon > ron; ARLov€ [ARLor, ARLoy |. 

Step-4: Generate a random sample of size ny«e [nr , ny] from the neutrosophic standard normal distribution 
with uine[p11, 41U] and variance ofe [o of]. Compute Xine |Xir, Xiu | for ith subgroup. 

Step-5: Compute the statistic MAjne[MAjz,MAjy] and plot it on LCLye[LCL¡, LCLy] and 
UCLye[UCL¡, UCLy]. Note the first out-of-control value, which is called the run length for the shifted process. 
Step-3: Repeat the process 10,000 times and compute the NARL and NSD, say ARL¡ye[ARLoz, ARLou] at 


Lin€[L1L> Miu] for various values of c. 
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0 [201.71, 211.83] | [199.03, 205.38] 


[297.65, 303.3] 


[289.36, 290.94] 


[369.7, 372.87] 


[349.46, 353.55] 

















0.1 [138.75, 106.01] | [137.11, 102.25] | [197.63, 140.91] | [196.55, 137.01] | [248.96, 173.35] | [241.98, 171.21] 

0.2 [68.64, 39.16] [66.07, 35.54] [93.49, 48.82] [91.93, 45.06] [114.03, 57.83] [111.2, 55.96] 
[35.38, 18.11] [33.46, 14.98] [45.87, 21.63] [43.65, 18.14] [54.14, 23.92] [51.64, 20.41] 
[19.48, 10.4] [17.59, 7.13] [24.34, 11.88] [22.35, 8.59] [27.89, 12.6] [25.69, 9.37] 
[12.12, 1.83] [10.18, 2.65] [14.42, 7.4] [12.51, 4.64] [16.03, 8.26] [14.25, 4.86] 

0.6 [8.12, 1] [6.21, 0] [9.31, 1] [7.33, 0] [10.15, 1] [8.24, 0] 

0.8 [1.98, 1] [2.22, 0] [5.07, 1] [3.04, 0] [5.37, 1] [3.37, 0] 








Table 1. The values of NARL when n e [4, 6] and w e [3, 5]. 


























0.1 [119.76, 89.24] | [119.08, 85.27] | [173.58, 117.97] | [170.18, 113.52] | [208.68, 143.88] | [206.34, 141.12] 
0.2 [49.92, 30.56] [48.38, 27.17] [65.68, 37.53] [64.44, 33.93] [79.91, 41.55] [78.47, 39.11] 
0.3 [23.42, 13.69] [21.36, 10.39] [29.52, 15.76] [27.77, 12.58] [34.03, 17.18] [31.75, 13.78] 
0.4 [12.56, 8.33] [10.68, 4.95] [15.18, 8.85] [13.17, 5.47] [16.89, 9.53] [14.72, 6.25] 
0.5 [7.84, 1] [5.87, 0] [9.05, 1] [7.04, 0] [9.77 1] [7.85, 0] 

0.6 [5.52, 1] [3.53, 0] [6.02, 1] [4.07, 0] [6.42, 1] [4.45, 0] 




















Table 2. The values of NARL when n e [6, 8] and we [3, 5]. 


[202.51, 211.04] 


[202.12, 204.53] 


[296.33, 300.31] 


[293.46, 291.44] 


[371.81, 371.09] 


[354.74, 353.76] 





[104.63, 78.17] 
[38.76, 24.34] 


[101.61, 74.31] 
[36.44, 20.94] 


[149.81, 102.88] 
[51.19, 29.23] 


[147.06, 100.28] 
[49.49, 25.84] 


[182.36, 122.74] 
[62.12, 32.76] 


[181.09, 118.33] 
[59.99, 29.49] 





[17.06, 11.25] 
[9.26, 1.01] 


[15.11, 8.06] 
[7.39, 0.2] 


[21.36, 12.63] 
[10.85, 3.39] 


[19.14, 9.33] 


[24.2, 13.62] 
[11.85, 7.05] 


[22.04, 10.26] 
[9.94, 4.57] 








0.5 [5.87, 1] 





[3.92, 0] 


[6.44, 1] 





[4.44, 0] 


[7.1, 1] 





[5.21, 0] 








Table 3. The values of NARL when n e [8, 10] and w e [3, 5]. 


[204.92, 208.42] 
[62.86, 42.28] 


[203.64, 204.08] 
[60.87, 38.62] 


[300.03, 295.65] 
[85.53, 52.85] 


[289.05, 286.77] 
[83.83, 49.27] 


[371.96, 375.79] 
[103.05, 61.52] 


[350.14, 353.69] 
[101.32, 58.7] 





[17.53, 11.27] 


[15.87, 8.1] 


[21.25, 13] 


[19.31, 9.79] 


[24.32, 13.93] 


[22.7, 10.69] 





[7.3, 1] 





[5.35, 0] 


[8.2, 1] 





[6.16, 0] 


Table 4. The values of NARL when n e [18, 22] and we [3, 5]. 
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[8.89, 1] 





[7.09, 0] 
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n= [4, 6] 
Existing chart Proposed chart 
C NARL NSD NARL NSD 
-= 1.5 [1, 1] [0, 0] [2.05, 1.33] [1.47, 0.67] 


[369.7, 372.87] | [349.46, 353.55] | [380, 370.71] [362.71, 353.45] 
[248.96, 173.35] | [241.98, 171.21] | [320.91, 294.43] | [313.19, 289.75] 
[114.03, 57.83] | [111.2, 55.96] [207.76, 164.71] | [207.45, 163.89] 
[54.14, 23.92] [51.64, 20.41] [126.96, 89.1 [127.24, 89.17] 
[27.89, 12.6] [25.69, 9.37] [73.79, 48.07 [73.14, 47.34] 
[16.03, 8.26] | [14.25, 4.86] [45.39, 26.98 [45.05, 25.99] 
[10.15, 1] [8.24, 0] [29.17, 16.53 [28.85, 15.95] 
0.8 [5.37, 1] | [3.37, 0] [12.75, 6.9] [12.26, 6.37] 


























Table 5. The comparison of the proposed chart and existing chart when NARL = 370. 


Using the NMC simulation process, ARLiye[ARL¡L, ARL¡y] and NSD for various c, nye[nr, nu] and 
WN€[wr, wu]are determined and placed in Tables 1, 2, 3 and 4. From Tables 1, 2, 3 and 4, the following trends 
can be noted in the values of ARL¡ye[ARL;¡, ARLiy]. 


1. For the fixed values of wwe [wr , wu | the values of ARL¡ye[ARL;¡1, ARL¡ y ] decreases as nye [nz, ny | increases. 
2. The values of ARL¡ye[ARL¡r, ARLiy ] increases as roy increases from 300 to 370. 


Results and discussion 

In this section, we will compare the performance of the proposed NMA chart over the existing MA chart under 
classical statistics in terms of ARL¡ye[ARLjz, ARL¡y]and by a simulation study. For this comparison, the same 
values of ron€ [ror rou |, €; wN€ [wr, wv] and nye[n,, ny ] are considered. 


Comparison in NARL. The proposed control chart is the extension of the X-bar control chart under neu- 
trosophic statistics proposed by Aslam and Khan”. The proposed chart reduces to Aslam and Khan” chart when 
WN€[1, 1]. We presented the values of ARL¡ye[ARL¡L, ARL1u]and NSD of both charts when roy e [370, 370] in 
Table 5. 

From Table 5, it can be noted that the proposed chart has smaller indeterminacy intervals of 
ARL\ne€[ARL 1, ARL¡y]and NSD as compared to the chart proposed by Aslam and Khan”. For example, when 
c=0.1, the values of ARL¡ye[ARL;¡1) ARL¡y]and NSD from the proposed charts are [248.96, 173.35] and [241.98, 
171.21], respectively. On the other hand, the values of ARL¡ye[ARL¡z, ARL¡y]and NSD from the existing control 
chart are [320.91, 294.43] and [313.19, 289.75], respectively. From this comparison, it is concluded that when 
c — 0.1, the proposed chart indicates the shift in the process between 173rd and 248th sample while the existing 
chart is expected to detect the shift between 294th and the 320th sample. We note that the proposed chart is an 
efficient chart in detecting the quick shift in the process as compared to Aslam and Khan” chart. 


Comparison by Simulation Data. Now we compare the performance of the proposed chart with the 
chart proposed by Aslam and Khan” and X-bar control chart under classical statistics using the simulated data. 
Let nye[6, 8], wy €[3, 5] and roye[370, 370] The 20 observations are generated from the neutrosophic normal 
distribution with uye[0, 0] and variance cie [1, 1]. Later 20 observations are generated from the process with 
a shift of 0.30. The values of statistic MA¡ye[MA¡L, MAjy ] and two existing charts are calculated and plotted in 
Fig. 1. From Fig. 1, it can be seen that the proposed chart detects shift at the 20th sample. On the other hand, the 
chart proposed by Aslam and Khan"? detects a shift in the process at around 30th sample while the traditional 
Shewhart control chart under classical statistics does not show any shift in the process. We also note that the 
proposed chart shows several points in the indeterminacy interval of control limits. 


Case study 

A meteorologist is interested to apply the proposed control chart for the better prediction and analysis of the 
weather. To explain the implementation of the proposed control chart, we collected the weather data in Jed- 
dah, Saudi Arabia. The temperature data of October 2019 at three different occasions is collected from https:// 
www.timeanddate.com/weather/saudi-arabia/jeddah/historic. The temperature data at morning 6, afternoon and 
evening 6 are reported in Table 6. From Table 6, it can be seen that the weather data has neutrosophic numbers. 
Therefore, the existing control chart under classical statistics cannot be applied for this type of data. The use of 
the existing chart under classical statistics may mislead the meteorologist. The proposed control chart under 
neutrosophic statistics can be applied for the monitoring of the weather. The proposed control chart and existing 
control charts for the weather data are shown in Fig. 2. 
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Table 6. The weather data. 
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Figure 2. The control charts for real data. Left figure (proposed chart) and right figure (Shewhart chart). This 
figure is generated using R version 3.2.1 (https://www.r-project.org/). 


From Fig. 2, we note that several temperature values are close to the upper value of the control limits. On the 
other hand, the existing chart under the classical statistic shows that the weather is in-control and the meteor- 
ologist needs no action. From this comparison, it is concluded that the proposed chart indicates some issue in 
weather and meteorologist should be alert about it. 
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Concluding remarks 

A new moving average control chart for the normal distribution was offered under the neutrosophic statistics. 
[he parameters of the offered chart were determined through simulation under neutrosophic statistics. The 
comparison study showed the superiority of the proposed chart over the moving average control chart under 
classical statistics. A real example from the weather was chosen to present the implementation of the chart. From 
the simulation study and real data, the proposed chart was found to be effective to be applied for temperature 
monitoring than the existing control chart. The proposed control chart can be applied in the weather department 
for the monitoring of the process. The proposed chart using EWMA statistics can be extended for future research. 


Data availability 


[he data is selected from https://www.timeanddate.com/weather/saudi-arabia/jeddah/historic. 
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